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Abstract 



The BUU transport model is applied to study strangeness and pion production 
in nucleus-nucleus collisions. Starting from proton induced reactions, we further 
investigate Si+Au, Au+Au and Pb+Pb collisions in the energy range between 
2 and 40 A-GeV and compare with data and with other transport calculations. 
The qq annihilation, or resonance, channel simulated by the string model in 
meson-nucleon collisions at y/s > 2 GeV is introduced. The importance of this 
channel for a good description of the proton-nucleus data on K + production is 
demonstrated. We, furthermore, show that the meson-meson collisions contribute 
significantly to the KK production in heavy-ion collisions above 5 A GeV and 
improve an agreement with data on the K + /tt + ratio. Finally, we study the 
influence of the in-medium modifications of the FRITIOF model on the pion and 
kaon production. 



1 Supported by GSI Darmstadt 

2 On leave from RRC "I.V. Kurchatov Institute", 123182 Moscow, Russia 



1 Introduction 



High-energy heavy-ion collisions offer a unique possibility to study nuclear matter 
at high densities and temperatures under laboratory conditions. The maximum 
compression is expected at the beam energy around 30 A-GeV. The most intrigu- 
ing phenomenon which can happen in highly compressed nuclear matter is the 
transition to the quark-gluon plasma. Pioneering work at large baryon densities 
was done at the AGS in Brookhaven where the energy range up to 15 A-GeV 
was explored (see [TJ and refs. therein). The future facility at GSI will provide 
beams from 2 A-GeV up to 40 A-GeV. One believes that an indirect signal for the 
quark-gluon plasma is the strangeness enhancement, which was first suggested 
by Rafelski and Miiller [2]. The enhancement should then be seen in the most 
abundant strange particles, the kaons. At the AGS and the SPS energies the 
K + /ir + ratio was studied and, indeed, a maximum in the ratio was found at 
about 30 A-GeV [B]. Theoretical calculations with different transport codes - 
RQMD jH 01, HSD |3] and UrQMD — have recently been performed. These 
calculations were not able to reproduce the peak in the ratio, which was either 
at a wrong energy (RQMD, UrQMD) or not present at all (HSD). In the case 
of HSD and UrQMD this discrepancy was due to overpredicted pion production, 
while the kaon yield was well described With an eye on the planned CBM 
experiment at GSI it is therefore important to check whether the mentioned dif- 
ficulties are genuine difficulties in the transport approach or consequences of a 
particular numerical implementation. 

In the present work we study pion and kaon production at the energies 2 - 
40 A-GeV within the BUU model (ZHSj; we stress that this is a numerical imple- 
mentation independent of those employed in jHlSEl- First, we systematically 
increase the system size and show its effect on particle production. In partic- 
ular, an analysis of the centrality dependence of the pion and kaon production 
from Au+Au collisions at 10.7 A-GeV is performed in comparison to the data 
Then we study the K + /tt + ratio in central Au+Au and Pb+Pb collisions as 
a function of the beam energy and compare our results to the experimental data 
and to other models. Special emphasis is put on the strangeness production and 
we will show in detail the most important production mechanisms at different 
energies. 

The structure of the paper is as follows. In Section El we describe our BUU 
model. In Section El we study the influence of the system size on the pion and 
kaon production. In Section 0] we show the excitation functions of pions, kaons, A 
and £ hyperons in central heavy- ion collisions. In Section El we discuss a medium 
modification of the FRITIOF string model and its influence on pion and kaon 
production. The summary and outlook are given in Section El 



1 



2 The BUU model 



Our calculations are based on the BUU model described in Refs. [HE]- In the 
high-energy range (y/s > 2 GeV) we adopt the treatment of Falter et al [Ej. Thus 
we will drop details concentrating only on the main ingredients and modifications. 

The model treats a nucleus-nucleus collision explicitly in time as a sequence 
of baryon-baryon, meson-baryon and meson-meson collisions. If not specified 
explicitly, the calculation is always done in the cascade mode, i.e. particles prop- 
agate freely between the two-body collisions. The baryon-baryon collisions at the 
invariant energy ^/i < 2.6 GeV are treated via a resonance scenario, whereas 
at y/s > 2.6 GeV a FRITIOF string model ^1] is applied. In the case of the 
meson-baryon collisions the resonance (FRITIOF) model is used at ^/s < (>)2 
GeV. In the most part of calculations we use an energy dependent strangeness 
suppression factor from Ref. [T2*] : 

f 0.3 for v/i > 20GeV 
7 = ttA = I 0-4 for v/i < 5GeV 

[0.433 - ^y/slGeVy 1 for 5 GeV < < 20 GeV 

(1) 

Sometimes we will also apply an energy independent strangeness suppression 
factor 7 = 0.3 which will be mentioned explicitly in the text. 

The most important modifications are the implementation of the strangeness 
production channels in meson-meson collisions and the possibility for a baryon 
and a meson to annihilate into a resonance with the invariant mass more than 
2 GeV, whose decay is simulated by the string model. These two modifications 
will be explained in detail below. In Section 0] we show how the meson-meson 
collisions and the effective resonance channel influence the K + /tt + ratio (see Fig. 

USD- 

2.1 Meson-meson reactions 

Our BUU model explicitly propagates tt, rj, p, a, u, (ft, K and the i^*-mesons. Ch- 
armed mesons are also included, but they are not relevant for the energies under 
consideration and thus we will not mention them further. At the beam energy up 
to 2 A GeV production of mesons heavier than pions is negligible and the only 
relevant meson-meson channel of strangeness production is tttt «-» KK, which 
was included in the BUU model earlier |Sj. At higher energies heavier mesons 
are produced more abundantly and, therefore, one also has to take into account 
the strangeness production in other meson-meson collisions, e.g. irp <-> KK. The 
problem is, however, that the cross sections of these processes are not measured 
experimentally. In Ref. the cross sections of the processes tttt — ► KK, np — ► 
KK and pp — ► KK were calculated. We use a parametrisation of nn — > KK 
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from Ref. which is based on these calculations: 

W = C 6.075 ( mb ), (2) 

where the factor C is the combination of the Clebsch-Gordan coefficients for the 
respective isospin channels: 

C = | (i\iim\mi \ i^IM) | 2 | (i 3 i 4 m 3 m4 | i 3 i 4 JM) | 2 , (3) 

7=0,1 

where ik and are the total isospin and the third isospin component of the 
particle k. Incoming and outgoing particles are enumerated by the pair of indices 
1,2 and 3,4 respectively. For simplicity we take the same cross section for pp — > 
KK and up — > KK since the isospins of the incoming particles are the same; 
this is not exactly the result from Ref. |13j . For all other reactions with two 
nonstrange mesons in the incoming channel we assume a constant value of 2 mb 
for the cross section. The back reactions are included and their cross sections are 
calculated according to detailed balance. 

By setting the cross sections constant we did not take into account the p-wave 
suppression of e.g. up — > KK near threshold. Due to the spin of the p-meson 
the strangeness production in that reaction is suppressed up to the threshold of 
up — > KK* and irp — > K K. Including this suppression, however, has only a 
small effect on our results (see discussions in section 4 and Fig l2*U|) . 

Since the elementary reactions are not accessible experimentally, we test our 
choice of the cross sections by looking on the yields of kaons in heavy-ion colli- 
sions. Although this will give us only a rough estimate of our cross sections due 
to presence of the baryon-baryon and the meson-baryon channels of the kaon pro- 
duction, there is no other way to get more reliable cross sections for meson-meson 
reactions. 

2.2 Annihilation processes 

Fig. d (see dashed lines) shows 3 that we underestimate the strangeness produc- 
tion in the region just above the FRITIOF threshold (a/s > 2 GeV). We see that 
directly above the threshold the cross section without the annihilation descends 
to almost zero. 

This is due to the fact that the FRITIOF model is only capable to produce two 
excited hadrons, which fragment separately. Thus it is not possible to describe 

3 In Figs. I1I3I5I7I the statistical errorbars are shown for the theoretical curves. They are 
calculated assuming a Poisson distribution, i.e. by dividing the plotted physical value by y/N, 
where N is the total accumulated number of events which is used to construct the value. In 
all other figures the statistical errors of theoretical results are either negligibly small or visible 
from the histogram representation. 
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Drell-Yan like processes in which, e.g., a quark from an incoming baryon and an 
antiquark from an incoming meson annihilate (c.f. Fig. |2J. An example for such 
a process is ttN — > YK, where Y stands for a hyperon. 

For that reason we have included the annihilation channel phenomenologically. 
In the case of a reaction of a baryon with a meson we check if an annihilation 
between a quark and an antiquark is possible; we split each of the interacting 
hadrons into their constituents and check whether a quark and antiquark with 
the same flavour exist. If they exist we annihilate the quark and antiquark 
with probability (jlj) and (jSJ, neglecting any particles that might be created in 
this process. In order to make up for this neglect we put all the energy and 
momentum of both incoming hadrons into the remaining quark content of the 
baryon and the meson. The fragmentation of this hadron is then done according 
to the Lund Model. UrQMD jT^j interprets meson-baryon reactions in a similar 
way. RQMD PB] also includes baryon resonances with mass > 2 GeV whose 
decay is described by the Lund model. 

The probability for the annihilation is chosen such that we agree with exper- 
iment for the strangeness production in irp — > strange particles (see solid lines in 
Fig. EJ): 

s/s 

Prob (annihilation) = max(0.85 - 0.17 • 0). (4) 

For the constant strangeness suppression factor 7 = 0.3 the probability for the 
annihilation processes is readjusted: 

fs 

Prob(annihilation) 7=0 3 = max(1.2 - 0.2 ■ ^— , 0). (5) 

There are two main reasons for the increase of the strange particle production by 
including the annihilation. First, we include new channels, as discussed above. 
Another point is that the invariant energy per string decay is higher. If we have 
two strings instead of one, the two strings decay separately and therefore it will 
often occur that every string alone is below the threshold for strangeness pro- 
duction. By putting all the energy into one string, the invariant energy becomes 
higher and the production of strangeness more probable. 



3 System-size dependence 

In order to clarify the reaction mechanisms in heavy-ion collisions, it is instruc- 
tive, first, to understand the proton and light ion induced reactions. For larger 
mass numbers of colliding nuclei, the effect of secondary hadron-hadron collisions 
becomes more and more important, which drives the system towards thermal 
equilibrium and enhances the maximum baryon density reached in the collision 
process. Thus by increasing the system size we can also see how particle spectra 
evolve with increasing density and if our model within the standard parameters is 
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able to reproduce the experimental measurements. In the discussion of numerical 
results it is assumed, if the opposite is not stated explicitly, that the meson-meson 
cross sections and the qq annihilation are included as described in the previous 
section. 

First, we study the proton- induced reactions p+Be and p+Au at the beam 
momentum of 14.6 GeV/c measured at BNL-AGS J7|. Fig. El shows rapidity 
distributions of produced tt^ and K ± for the p+Be system. In this case of light 
target, the pions and kaons are produced mostly in the first-chance nucleon- 
nucleon (NN) collisions and have only a small probability to rescatter afterwards. 
Therefore, their rapidity distributions are centered near the NN center-of-mass 
(cm.) rapidity yNN = Vbeam/^ = 1.72. The pion yield is underestimated by 
~ 20%, whereas the kaon and antikaon yields are well described by BUU for the 
p+Be system. 

According to Ref. [Ej, we have fitted the calculated transverse mass spectra 
with an exponential function 

d 2 a 

— = a exp{— m±/T} (6) 

27im±dmj_dy 

at various rapidities y. Fig. HJpresents the inverse slope parameter T as a function 
of rapidity for tt + and K + in the case of p+Be collisions. The pion inverse slope 
parameter is well reproduced by BUU except for the very forward and the very 
backward rapidities in the NN cm. system. The calculated kaon inverse slope 
parameter overestimates the data by ~ 20% at midrapidity. 

Fig. El shows the rapidity distributions of n ± and K ± for the p+Au collisions. 
These distributions are shifted to smaller rapidities y < yNN with respect to the 
case of p+Be reaction (Fig. El) due to contribution of the secondary NN and 
7rN collisions to the meson production and rescattering of the produced mesons 
on the target nucleons. The K~ rapidity distribution is narrower and is shifted 
somewhat less than the K + distribution, since an antikaon is always produced 
together with a kaon, while a kaon can be also produced in association with a 
hyperon which requires less cm. energy Thus, the secondary NN and 7rN 
collisions contribute more to the K + than to the K~ production. BUU describes 
the experimental pion and kaon rapidity distributions within ~ 20%. In Fig. E]we 
present the rapidity dependence of the inverse slope parameter T of the K + and 
7r + transverse mass spectra for the p+Au reaction. There is a good agreement 
between BUU and the data except for the very forward rapidity in the 7r + case 
where we overpredict the experiment by ~ 25%. 

In agreement with the data, we observe little change in the value (~ 150 MeV) 
of the inverse slope parameter for 7r + and K + with increasing target mass (cf. 
Figs. Eland El). The K + yield at y = yNN shows a factor of two enhancement, 
while the K~ yield at y = yNN stays practically unchanged both in BUU and 
in the data (cf. Figs. ElandEJ)- We attribute this behaviour to the stronger 
absorption of K~ in the heavier target. The experimental pion yield at y — yNN 
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is the same for both systems, whereas in BUU we observe a slight enhancement 
of the pion yield at y = i/nn with increasing target mass. 

Next, we present results for Si+Au collisions at the beam momentum of 14.6 
A GeV/c, which were studied experimentally in Ref. pQ. This reaction has earlier 
been studied theoretically in Refs. [TT)| H£] . In Ref. ^E] the RQMD model has 
been employed, which includes known nonstrange baryon resonances with m < 2 
GeV. RQMD also describes the production of baryon resonances with m > 2 
GeV in high-energy collisions; their decay is simulated by a string model. Shapes 
of the 7T and K transverse momentum spectra are quite well described by 
RQMD. No conclusions on the agreement of the absolute yields of the produced 
particles with data have been drawn in JH] due to the absolute normalization 
on the experimental 7r + spectra. However, ratios K + /tt + and K~/tt~ computed 
within RQMD agree well with data. In Ref. fH] a relativistic hadronic cascade 
(ARC) model has been used; a pure hadronic scenario without string excitation 
was assumed. A good agreement between ARC and E-802 data on proton m± 
spectra, proton, tt + and K + rapidity distributions was reached within the reso- 
nance model, i.e. when, e.g. a three-pion production channel in a NN collision is 
simulated as NN — > A An rather than directly as NN — > NNrr-k. Inverse slope 
parameters of the m± spectra for protons and pions are well described by ARC, 
but underestimated by ~ 20% for kaons. 

We have considered only central collisions Si+Au corresponding to 7% of the 
inelastic cross section selected on multiplicity of charged particles [T]. In the 
theoretical calculations we selected the central collisions in the same way. Fig. 
shows the calculated 7r =l= and K rapidity distributions which were divided by 
the projectile mass (28) in order to be able to directly compare them with the 
rapidity distributions from the proton induced reactions. 

By comparing the data points in Figs. [7| and 03 we see that the pion yields at 
V = Dnn are, practically, the same in p+Au and Si+Au systems. This feature is 
not reproduced by BUU: there is an enhancement of the pion yield per projectile 
nucleon in the system Si+Au with respect to the p+Au system in our calculations. 
This may indicate a problem with pion production (or reabsorption) in the heavy 
system. The experimental K yields per nucleon are higher in the Si+Au case 
than in the p+Au case, which is well reproduced by BUU. 

In Fig. |H]we present the inverse slope parameters of the K + and ir + transverse 
mass spectra. Despite of the big errorbars plus systematic errors of ±10 % which 
are not included into the errorbars of the experimental data |XJ|, we see that 
BUU underpredicts the inverse slope parameter for K +, s by about 25 % and for 
7r +, s by 15 %. The calculated inverse slope parameter stays, practically, constant 
T ~ 150 — 160 MeV for all three systems p+Be, p+Au and Si+Au both for 7r + 's 
and K +, s, whereas the experimental data show higher T ~ 200 MeV for K +, s in 
the Si+Au system. 

Studying strangeness production in more detail, we have also performed the 
BUU calculations using the constant energy-independent strangeness suppression 
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factor 7 = 0.3 (dashed lines in Figs. EIEIIZI)- The K + rapidity distributions favour 
the energy- dependent strangeness suppression factor, while the K~ spectra are 
better described with 7 = 0.3. 

In order to demonstrate an effect of the qq annihilation on the K + production 
(see discussion in the previous section) we also show in Fig. El the results without 
the annihilation. In the p+Au system the secondary irN and pN collisions play 
already an important role. Thus, including the annihilation improves an agree- 
ment with the data (c.f. solid and dotted lines in the lower left panel of Fig. 



The heaviest colliding system measured at AGS is Au+Au at the beam ener- 
gies of 2-10.7 A GeV [5J EH 1201- Before discussing the beam energy dependence 
(see next section), we will consider the centrality dependence of the pion and 
kaon production for Au+Au collisions at the top AGS energy of 10.7 A GeV 0. 

In Ref. [0] the collision centrality was determined by using two criteria: (i) The 
energy deposited in the zero-degree calorimeter Ezcal, which gives an estimate 
of the projectile participant number N pp as 



where Ej™ m = 2123 GeV is the kinetic energy of the beam. The smaller Ezcal 
is, the larger is the size of the participant zone, which selects geometrically more 
central events, (ii) The multiplicity of particles with velocity (3 > 0.8 in the New 
Multiplicity Array (NMA) mult N MA- The velocity cut filters out the slow pro- 
tons, whereas the produced mesons (mostly pions) are accepted. Thus, the larger 
multNMA corresponds to the larger energy transfer from the longitudinal motion 
of colliding nuclei to the meson production. In average, events with smaller im- 
pact parameter b have larger multNMA- However, at fixed b the multiplicity 
multNMA fluctuates stochastically event-by-event depending on the amount of 
stopping of the counterstreaming nuclear matter in the interaction zone. 

According to Ref. [H] we, first, divided BUU events using Ezcal into eight 
classes with increasing Ezcal (decreasing centrality) from the first to the eighth 
class (see Table II in 0). Second, by modelling the NMA acceptance, we sub- 
divided each of the first three E Z cal event classes to the three multNMA classes 
with decreasing multNMA (decreasing centrality) from the first to the third mul- 
tiplicity class (see Table III in 

Fig. EH shows the proton rapidity distributions for various combinations of 
Ezcal- and multNMA- event classes. We see that in each case BUU overestimates 
stopping. Agreement with experiment can be improved by taking into account 
the in-medium corrections to the FRITIOF model (see Sect V). 

Figs. ITUl and HP show the rapidity spectra of pions and kaons for different event 
classes selected by applying the Ezcal cut only. The K + rapidity spectra are very 
well described for all E Z cal classes. For the 7r + rapidity spectra we see deviations 




(7) 
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from the data: In the most central collisions there are excessive pions in BUU 
produced mainly at midrapidity. With decreasing centrality the deviation from 
the data disappears gradually, and in the most peripheral collisions there is even 
an underprediction of the pion multiplicity by BUU. These results are consistent 
with Figs. El El 13 where one can also observe a tendency to overpredict pion 
production with increasing size of the participant zone. 

Figs. ^1 and show transverse mass spectra of 7r + 's and K + 's for the 
central collisions of Au+Au at 10.7 A GeV. The spectra are shown for various 
rapidities starting from the backward rapidity in the cm. frame (upper line) to 
the midrapidity (lower line). The shapes of the 7r + spectra are well described by 
standard BUU, however, the 7r + yields are slightly overpredicted at small m± (see 
also upper left panel in Fig. HUjl . The agreement of BUU with the K + spectra 
is much poorer. The low m_|_ part of the K + spectra is, typically, overestimated 
by BUU, whereas the high m± part of the spectra is underestimated by our 
calculations. Thus, BUU underestimates the inverse slope parameter of the K + 
transverse mass spectra, while the K + yield is well described (c.f. upper left 
panel in Fig. fTTj) . This problem has been pointed out earlier in Ref. 123] . 

Fig. IT41 shows a fiducial yield of K + and tt + divided by the projectile partici- 
pant number N pp as a function of N pp . The fiducial yield is defined as follows 9J: 



where the dN/dy are the rapidity distributions selected by the zero degree energy. 
The K + fiducial yield, as expected, agrees quite well with the data except for 
a slight underprediction at peripheral collisions. The tt + fiducial yield increases 
with N pp faster than the data do. In the absence of secondary NN collisions the 
fiducial yields divided by N pp would be constant. 

4 Excitation functions 

In this section we show the excitation functions of pions, kaons, A- and E- 
hyperons from central Au+Au and Pb+Pb collisions in comparison to data, two 
other transport models HSD and UrQMD [H| and the statistical model [22| • The 
calculations performed with the transport models have all been done in the cas- 
cade mode, which makes the comparison easier. We selected the data sets for 
the Au+Au system at 1.96, 4.00, 5.93, 7.94 and 10.7 A GeV [131201 with 5 % of 
the most central events and for the Pb+Pb system at 30 and 40 A GeV [HI I2B] 
with 7 % of the most central events. In the theoretical calculations we used a 
sharp impact parameter cut off at 3.5 fm for AGS energies and 4 fm for SPS 
energies. The influence of the centrality selection was tested at 10.7 A GeV by 
comparing calculations with a sharp cutoff to calculations which were done by 





8 



employing the centrality criteria described in sectional No deviations were seen 
in the observables which will be presented in the following. 

Fig. "l"5l shows the midrapidity yield of positive pions as a function of the beam 
energy. We see that all three models overpredict the pion yield in the considered 
beam energy range (2-40 A GeV). Our model (solid line) overpredicts the n + 
midrapidity yield by ~ 10 % at 40 A GeV to ~ 50 % at 2 A GeV. However, the 
shape of the experimental excitation function dN / dy vs i?Lab fc> r 7r+ is remarkably 
well described by BUU. The HSD model (dot-dashed line) produces the 7r + yields 
close to the BUU results excepting the points at 6, 8 and 10.7 A GeV, where HSD 
has ~ 10 % more pions than BUU. The UrQMD model (dotted line) agrees well 
with the pion data at the smallest energy of 2 A GeV, but the pion yield grows 
too fast with energy within UrQMD producing the discrepancy ~ 30 % with data 
at the highest considered energy of 40 A GeV. 

Figs. [El IT7I and fTHl show the midrapidity yields of K + ,K~ and (A + E°), 
respectively, as functions of the beam energy. BUU quite well describes the K + 
midrapidity yield excepting the points at 4 and 6 A GeV, where BUU overesti- 
mates the data by 30-50 %. 

The K~ midrapidity yield and (A + E°) midrapidity yield at E Lah < 40 A GeV 
are overestimated by BUU. Using the constant strangeness suppression factor 
7 = 0.3 (dashed lines) reduces the yields of K + , K~ and hyperons. This leads to 
a better description of the K~ yields, while in the cases of K + and (A + £°) it is 
hard to judge which strangeness suppression factor works better. In Fig"i"""l and 
FigEHwe see that up to about 10 A GeV the K + and hyperon yields are better 
described with 7 = 0.3. Above that energy, however, the energy dependent 
suppression factor works better. For the antikaons the strangeness exchange 
processes KN <-> irY are important due to the strong in-medium modifications 
[2*4] . which are not taken into account in our study. 

The HSD results on the K + production are close to our calculation with 
7 = 0.3. The UrQMD model gives the K + yield at lower energies similar to our 
standard calculation, whereas at higher energies UrQMD produces substantially 
less kaons. The K~ yield is rather well described by both models, HSD and 
UrQMD, except for the point at 30 A GeV. At lower energies the (A + £°) yields 
calculated within HSD and UrQMD are somewhat closer to the data than our 
standard calculation. At higher energies the HSD, UrQMD and our standard 
calculation give very close results for the (A + E°) yield. 

Fig. EH shows the ratio of midrapidity yields of K + and 7r + as a function of 
the beam energy. In the upper left panel of Fig. flQlwe see that neither BUU nor 
HSD and UrQMD describe the ratio K + /tt + in the whole beam energy region. 
At the lowest beam energy of 2 A GeV BUU and UrQMD overpredict the ratio 
by a factor of two, whereas HSD agrees with data. Between 4 and 8 A GeV BUU 
is quite close to the data. However, the K + /tt + ratio excitation function levels off 
too early in BUU and, as a consequence, we underestimate the ratio by ~ 25 % 
between 10 and 30 A GeV. The HSD results on the ratio K + /tt + have a similar 
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beam energy dependence, but the value of the ratio is smaller by ~ 20 %, which 
is close to our calculation with 7 = 0.3. UrQMD produces a larger slope at lower 
energies overestimating the ratio at i^Lab < 8 A GeV, but at higher energies the 
slope gets negative which causes a strong discrepancy with data in SPS energy 
regime. Overall, we observe that BUU has the best agreement with data on the 
K + /tt + ratio in the considered energy regime. At the beam energies of 4-6 A 
GeV this, however, comes about due to cancellation of the overestimation of the 
pion and kaon yields. 

Fig. ED shows the ratio of the midrapidity yields K~ / K + vs the beam energy. 
The BUU calculations with the energy dependent strangeness suppression factor 
and our "standard" meson-meson cross sections is shown as the upper boundary 
of the errorband. In order to estimate an effect of the p-wave suppression on the 
up — > KK and ttuj — > KK cross section (see section |2~T|) . we also performed a 
calculation by putting the cross section equal to zero below the K*K production 
threshhold, which is shown by the lower boundary of the errorband in Fig]2"Ul 
Thus the p-wave suppression could reduce the K~ multiplicity by about 10%. 
The K + multiplicity is reduced less than 5% by this effect, since the fraction of 
K + coming from meson-meson collisions is less than the fraction of K~ coming 
from meson-meson reactions. (For this reason our results on the K + /tt + ratio 
FigllHlare, practically, untouched by the p-wave suppression effect.) 

Overall, we see that BUU overestimates the K~/K + ratio independent on 
the strangeness suppression factor. This result is expected from the previous 
Figs. and El where we see that the K + yield is rather well described by 
BUU, while the K~ yield is overestimated substantially. Since the choice 7 = 0.3 
reduces both K + and K~ yields, the ratio K~ / K + is practically independent on 
the strangeness suppression factor. The HSD and UrQMD models describe the 
experimental K~/K + ratio quite well. This can be also traced back to Figs. [TBI 
and EH 

Fig. I2T1 shows the ratio of the midrapidity yields (A + S°)/7r vs the beam 
energy. This ratio has a peak near i?Lab = 8 A GeV which is reproduced by BUU 
and HSD models. As far as the absolute values of this ratio are concerned at 
small energies E Lah < 10 A GeV, we would like to remind that both pion and 
hyperon yields are overestimated by standard BUU (see Figs. ITol and 1TS|) . Thus, 
the agreement of standard BUU with data on (A + S°)/7r at 8 and 10.7 A GeV 
is again a result of a mutual cancellation of the (A + S°) and ir excesses. The 
choice of 7 = 0.3 which describes the (A + S°) midrapidity yield at small energies 
better (c.f. Fig. 1T%|) leads to the underestimation of the (A + T l °)/n ratio at the 
peak due to the overestimation of the pion yield. 

All calculations discussed above were performed in the cascade mode. There 
is an option in our BUU model, which switches on a nuclear mean field potential. 
The nuclear mean field potential is necessary, in particular, for a description of 
the experimental data on collective in-plane and out-of-plane proton and neutron 
flows at 0.15-2 A GeV. At higher energies, however, the parametrisation of 



10 



the momentum-dependent interaction used in "21"" leads to too strong repulsive 
in-plane flow (see also Ref. [26 ). Nevertheless, in order to estimate the mean field 
effect on pion and kaon production, we have also done the calculation with the 
mean field potential (incompressibility K = 215 MeV, soft momentum-dependent 
mean field SM). The lower left panel of Fig. fTHl shows the results of this calculation 
(dashed line) in comparison with our standard calculation in the cascade mode 
(solid line) and with experimental data. We see that at E Lab < 40 A GeV the ratio 
K + 1 7r + is reduced due to the mean field potential, since the pion yield is relatively 
insensitive to the nuclear mean field, whereas the kaon yield is reduced. Indeed, a 
part of the kinetic energy of the counterstreaming nucleon flows transforms now 
to the potential energy. This reduces, generally, particle production. However, 
the kaon production is closer to its threshold than the pion production. Therefore, 
kaons are more strongly influenced by the mean field potential than pions. 

As we described in Section [21 the meson-meson collisions and the qq annihila- 
tion channel for the meson-baryon collisions are implemented in our BUU model. 
The lower right panel of Fig. []"""] shows an effect of these implementations on 
the K + /tc + midrapidity ratio. The result of our standard calculation including 
both the meson-meson collisions and the qq annihilation is shown by the solid 
line in Fig. [""""J The dotted and dashed lines represent the calculations without 
the meson-meson collisions but with annihilation and without the annihilation 
but with the meson-meson collisions, respectively. The meson-meson collisions 
strongly enhance the K + /tt + ratio above 6 A GeV due to the increased KK 
production. An effect of the qq annihilation channel is less pronounced: only a 
slight enhancement of the K + /n + ratio is visible at 5-10 A GeV. 

The upper right panel of Fig. El compares the BUU and the statistical model 
[2*2*] calculations for the K + /rr + ratio at midrapidity. Since we use a string model, 
which produces a multiparticle final state for the two colliding particles, the 
thermal equilibrium would be only reached if we also included the corresponding 
back reactions (c.f. Ref. [22]). However, this was out of scope of the present 
work. Nevertheless, there is a surprisingly good agreement between BUU and 
the statistical model. The statistical model is closer to the data, but it is also 
unable to describe the data points at 10 and 30 A GeV. 

Finally, in Fig. [221 we show the inverse slope parameter T of the K + transverse 
mass spectra vs the laboratory energy. At the energies 2, 4, 6 and 8 A GeV the 
inverse slope parameter was obtained by fitting the exponential function (j""J) to the 
m_i_-spectrum of kaons in the rapidity range \ (y — Dnn) /vnn\ < 0.25 20J. At 10.7 
A GeV the rapidity range was \(y — Vnn) Ivnn\ < 0.125 [20]. At the SPS energies 
of 30 and 40 A GeV the rapidity range was taken as \y — yNN\ < 0.1 [3I2S]- Our 
model underestimates the inverse slope parameter by 30-40 % (see also Fig. fT3*|) . 
A similar problem was reported before in Ref. ~2~T| "28] and ascribed to the lack of 
pressure due to missed nonhadronic degrees of freedom in the transport models. 
We speculate here that also the inclusion of multi-baryon collisions would tend to 
make the spectrum harder. The probability for such processes naturally increases 
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with a high power of baryon density. 

In order to see an origin of the produced kaons we performed a channel de- 
composition of the s-quark production for the central Au+Au collisions at 4, 10.7 
and 20 A GeV. Fig. ESI shows the number of the produced s-quarks vs time for 
four different channels: (i) baryon-baryon channel at high energies (y/s > 2.6 
GeV) or, in other words, the baryon-baryon reactions simulated by the FRITIOF 
string model (solid line), (ii) baryon- meson collisions at \fs > 2 GeV simulated 
by the string model (dashed line), (iii) baryon-meson collisions at y/s < 2 GeV, 
i.e. below the string model threshold (dash-dotted line), (iv) meson-meson colli- 
sions (dotted line). The low energy (y/s < 2.6 GeV) baryon-baryon collisions do 
not contribute to the s-production significantly at the considered beam energies. 
Thus, this channel is not shown in Fig. |2H1 We counted only the creation of the 
s-quark and we did not consider reactions or decays, as e.g. K* — > Kn where the 
s-quark is only shifted from a K* to a K. 

The baryon-baryon-string channel plays the dominant role in the whole beam 
energy region. This channel includes mainly the first-chance NN collisions be- 
tween the projectile and the target nucleons. The meson-meson channel is not 
important at 4 A GeV, but its role grows quickly with energy and at 20 A GeV 
it includes already 25 % of the produced s-quarks. A relative contribution of 
the baryon-meson-string channel also increases with energy, while the low energy 
baryon-meson collision relative contribution stays always very small and decreases 
with energy. 

The time evolution of the s-quark production can be better understood if one 
looks also on the central density time evolution shown in Fig. E]for the central 
Au+Au collision at 10.7 A GeV. The central density reaches its maximum value 
~ 4.5po &t t — 10.5 fm/c and stays above 3po i n the time interval t = 7.5 — 14 
fm/c, where the s-quark production just takes place. Here po = 0.16 fm~ 3 is the 
nuclear saturation density. Thus, strangeness is produced during the high-density 
stage of a heavy-ion collision. It is evident also from Fig. |2Hl that the s-quark 
production from the meson-meson and the baryon-meson channels, which contain 
the secondary collisions, starts later than from the baryon-baryon channel. 

5 In-medium modification of the FRITIOF 
model 

In the course of a heavy-ion collision the elementary hadron-hadron collisions 
happen at a finite baryon density. Therefore, the wave functions of incoming and 
outgoing particles are the in-medium plane waves rather than the vacuum ones. 4 

4 Exchange particles expressed by the propagators get also in-medium modified. However, 
this last effect would be strongly dependent on the model used for the description of an ele- 
mentary collision. Thus we will, for simplicity neglect it by using the vacuum matrix elements 
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In order to take into account the in-medium modifications of the incoming 
and outgoing particles in the FRITIOF events, we follow here the approach of 29 
generalized to the processes with many-meson final states. Only the events with 
two colliding nonstrange baryons will be modified. The in-medium modifications 
of the meson-baryon and meson- meson collisions are neglected, since as we expect, 
they are small with respect to the baryon-baryon case (see below). 

Let us consider the process 

B 1 B 2 -> B 3 B 4 M 5 M 6 ...M N , (9) 

where Bi,B 2 and -63,-64 are the incoming and the outgoing baryons, respectively; 
M 5 ,M 6 ,...,Mn are the produced mesons. The in- medium differential cross section 
of this process is given by the following expression 

d(jmed = M 4 (2m;)(2m;K2m5)(2m;) 1 ^ d ^_ 2(p; + ^ ^ ^ ^ ^ 

__ (10) 
where |T| 2 is the matrix element squared in the normalization of Ref. [3U] av- 
eraged over spins of initial particles and summed over spins of final particles, 

d^ N - 2 {pt + Pa! pIpI k* 5 , k*, k* N ) = 5 (4 \pl +p*-p* 3 -pl-kt-k*-...- k* N ) 

d 3 p* d 3 pl d 3 k* d 3 k* N 

X (W (2vr) 3 2(^)o (27r)32(fc 5 *)o ' ' ' (2ir)*2(k' N )° 

(11) 

is the (iV-2)-body phase space gl] with (p*)° = ((p*) 2 + (m*) 2 ) 1 / 2 , i= 1,2,3,4 
and (k*)° = ((k*) 2 + (m*) 2 ) 1 / 2 , i = 5, 6, ...,N being the zeroth components of 
the kinetic four-momenta, and the m* being the effective masses of the particles 
involved. In Eq. (|10p the flux factor is 

I* = q(y/^*, m*, m* 2 )Vs* , (12) 

where s* = (p* + p* 2 ) 2 and 

gf^K,™;) = [(5* + K) 2 - K) 2 ) 2 /(4 S *) - K) 2 ] 1 / 2 (13) 

is the cm. momentum of incoming baryons. 

The matrix element |T| 2 entering into Eq. (p[T)|) can be extracted from the 
vacuum cross section by dividing out the vacuum phase space and multiplying 
by the vacuum flux factor. Thus, our final result for the in-medium total cross 
section of the process is 

a med {V^) = Fa vac {^) . (14) 
on the place of the in-medium ones. 
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The modification factor F is 



F _ m\m\m\ml I <5> N _ 2 (^/¥; m*, m*, m* N ) 
m 1 m 2 m 3 m i I* <3>at_ 2 (\/s; rn 3 , m 4 , m N ) 

where / = q(y/~s,mi,m2)\fs- In Eq. (|T4j) yfs is the cm. energy of the colliding 
baryons in vacuum, which is directly provided by the BUU calculations in the 
cascade modus performed in the present work. The in-medium cm. energy \rs* 
is then determined from the condition that the energy excess above threshold is 
the same as in vacuum, i.e. 

VS* = a/s — m 3 — m 4 — ■ • • — m N + m*^ + ml + ■ ■ ■ + m* N . (16) 

Since the modification factor F is proportional to the product of the ratios of the 
Dirac mass to the bare mass for incoming and outgoing fermions, we expect that 
the meson-baryon and meson-meson channels will be modified relatively weaker. 

In Eq. (jl()j) we replaced the canonical four- momenta by the kinetic ones in the 
^-function entering the phase-space volume element ([lip. This is possible only 
if the vector fields cancel each other, which is valid in the case of the Walecka 
model (c.f. ref.[32j), but would be violated in a more sophisticated relativistic 
mean field model with momentum-dependent scalar and vector fields [33J. Taking 
into account the momentum dependence, in particular, for the vector field, which 
drops with momentum, is important for the description of the baryon flow in 
heavy ion collisions above 1 A GeV [31]. However, in the present exploratory 
work we will neglect the momentum dependence of the a and u fields, which 
would strongly complicate the calculation of the phase space volume. 

We evaluate the in-medium masses using a nonlinear version NL2 [HE] of the 
relativistic mean field model and assuming that the nucleons and all nonstrange 
baryonic resonances are coupled to the scalar mean field a and to the vector mean 
field us by the same universal coupling constants g a and g u |36^. This gives the 
Dirac effective masses 

m* B =m B + g a a (17) 

and the kinetic four-momenta 



P*b=Pb~ QluUJ (18) 

of the nonstrange baryons. The hyperon coupling constants are (c.f. Ref. |32j ) : 

2 2 

9a = 3^ , gl = g&>, (19) 

where Y = A or S. The baryon single-particle energy is 

e(pfl) = W + yJ(p* B ) 2 + (m* B ) 2 . (20) 
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For the mesons tt, p and u we neglect any in-medium modifications, while for the 
K and K single-particle energies we use the model of Ref. jSZj with parameters 
ofRef. [HE]: 

u(k,p) = ^/(k*y + (m* K y±V°, (21) 
where the upper (lower) sign corresponds to the K (K) case, 

k* = k =f V (22) 

is the kaon kinetic momentum, 



m*K= \h 2 K-^jfPs + V* (23) 



is the kaon effective (Dirac) mass, and 



Kf: 

is the kaon vector field. p s and j M are the scalar density and the baryonic four- 
current, respectively The parameters which appear in Eqs. fl2HI),(J23i) are T,kn = 
450 MeV, f n = 93 MeV and (f*) 2 = 0.6(U) 2 M- Within these parameters the 
following relation is expected to hold 



V* ~ l -g„u» . (25) 

Taking into account relations and (|25jl one can see that the vector field is, 
indeed, completely excluded from the energy-momentum conservation conditions 
for the strange particle production processes like B\B 2 — > B 3 Y±K or BiB 2 — > 
B^B^KK, which gives a possibility to simplify the in- medium calculations by 
just replacing the bare masses of particles by the Dirac masses and canonical 
four-momenta by the kinetic four-momenta. 

We have calculated the modification factor F(y/s,p) as a function of the 
cm. energy y/s and the baryon density p for various outgoing channels with 
no more than four mesons in the final state. We assume that the incoming 
baryons are nucleons, but an outgoing baryon can be either nucleon or A-hyperon. 
For outgoing mesons we have considered all possible combinations of pions, p- 
mesons, kaons and antikaons with no more than one kaon and one antikaon in 
the final state. The upper panel of Fig. 12*51 shows the medium modification factor 
F(y/s,po) for some selected processes: NN -> iViVTr, NNtttc, NNp, NAK, 
NNKK. We see that the modification factor depends on the outgoing channel 
rather weakly. In particular, the addition of a pion does not change the factor. 
On the lower panel of Fig. EBJ we demonstrate the density dependence of the 
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modification factor for the one-pion production channel. One can observe a strong 
decrease of the factor with the baryon density. 

For the application to the FRITIOF model built in the BUU code, we have 
stored the factors F on a two-dimensional grid (y/s,p). Once some final state is 
generated by FRITIOF, it is accepted with the probability F. In the case where 
at least one of incoming baryons is a resonance, we use the modification factor for 
incoming nucleons at the same \fs. If an outgoing baryon is the A-resonance, the 
factor for an ougoing nucleon is applied, shifted by the production threshold, i.e. 
F(y/s — ira + mjv). An analogous threshold correction is performed if S-hyperon 
and/or K* are produced. 

In order to see an effect of the in-medium string model modifications in heavy- 
ion collisions, we have performed the calculation for the Au+Au system at 10.7 
A GeV. The results of this calculation are shown by the dashed lines in Figs. 

The proton rapidity spectra (Fig. EJ) get now wider, in agreement with the 
experiment. This is expected since the inelastic NN cross section is reduced by 
the in-medium effects. The results with in-medium modifications for the most 
central events (upper left panel in Fig. EJ) have a big statistical error due to a 
small number of events in this centrality class. 

The 7r + rapidity spectra (Fig. ITUj) are reduced in a closer agreement with the 
data, except for the very peripheral collisions. The K + rapidity spectra (Fig. ITTf 
are also reduced. This, however, makes the agreement with the K + data worse. 

The transverse mass spectra of 7r + 's (Fig. IT2j) and K +, s (Fig. IT3|) are reduced 
more at high transverse masses, since the high-mj_ particles are emitted from 
hard collisions which happen at the high baryon density where the in-medium 
modifications are stronger. Thus, the m^-spectra get steeper, which, again, leads 
to a worse description of the K + data. The 7r + ' transverse mass spectra at small 
m_i_ are now better described, whereas at high m_|_ the in-medium modifications 
result in a slight underestimation of the experiment. 

Fig. E] summarizes our findings on the in-medium modifications. The ir + 
fiducial yield is described better with the in-medium modifications, except for 
the very peripheral collisions. The K + fiducial yield is underestimated at all 
collision centralities with the in-medium modifications. 

6 Summary and outlook 

In the present work, the BUU model developed earlier in Refs. [7J |H] is further 
improved by including the heavy (m > 2 GeV) resonance or gg-annihilation 
channel in the meson-baryon collisions and by the new meson-meson channels 
(np — > KK, pp — > KK) for the strangeness production. Moreover, an in-medium 
modification of the FRITIOF model by taking into account the effective (Dirac) 
masses of the incoming and outgoing particles is formulated and implemented in 
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BUU. The BUU model is applied to the nucleus-nucleus collisions at 2-40 A GeV. 

By performing the systematic study of the pion and kaon production for 
various systems and collision energies we came to the following conclusions: 

1. The and rapidity spectra and the inverse slope parameters of the 
transverse mass spectra from the proton-nucleus reactions are well described 
by BUU. The gg-annihilation channel improves the agreement with the data 
on the K + production in proton- nucleus collisions (Fig. EJ. 

2. In the peripheral Au+Au collisions at 10.7 A GeV the 7r + rapidity spec- 
tra are slightly underestimated, whereas the K + rapidity spectra are well 
described. This result is consistent with a good agreement between BUU 
and the data on the proton-nucleus reactions. The proton rapidity spectra 
are too narrow, i.e. the stopping power is overestimated by BUU, for all 
collision centralities. 

3. In the central nucleus-nucleus collisions (Si+Au at 14.6 A GeV/c, Au+Au 
at 2-10.7 A GeV, Pb+Pb at 30 and 40 A GeV) the K+ yields are, overall, 
well described, the K~ and hyperon yields are somewhat overestimated, 
the 7r + yields are overestimated in all central heavy-ion collisions under 
study. The inverse slope parameters of the K + transverse mass spectra are 
strongly underestimated (Figs. IT3*1 l2*2* j) . The pion slopes are well described 
(Fig. H2J). Our BUU results on n + and K + agree, generally, with results 
of the HSD and UrQMD calculations from Refs. [HI E] . The excitation 
function of the K + /ir + ratio is described, however, better and closer to 
the thermal model results due to the introduction of the new meson-meson 
channels of the strangeness production which increase the K + yield at the 
beam energy above 6 A GeV (Fig. ITT?j) . 

4. The in-medium modification of the FRITIOF model, which reduces the 
particle production cross sections, was tested for the system Au+Au at 
10.7 A GeV. This leads to the better description of the pion production, 
whereas the K + production is overdamped by the in-medium effects. The 
stopping power of the nuclear matter is reduced, which results in the better 
description of the proton rapidity spectra. 

The K + /tt + ratio is quite sensitive to the meson- meson cross sections and, 
within the reasonable choice of these cross sections, we have decreased the dis- 
crepancy between BUU and the data by about a factor of 2. The resulting 
disagreement ~ 20 % can, thus, hardly to be considered cLS cL SI eiial of the "new 
physics". The inverse slope parameter of the K + transverse mass spectra is a 
more serious problem for BUU. It would be worthwhile to study this topic in 
more detail by introducing, e.g. string-string and many-body collisions, which 
could both enhance the hard part of the kaon spectra. 
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Figure 1: The cross section of the strangeness production in 7r + p collisions (upper 
panel) and n~p collisions (lower panel) as a function of beam energy in compar- 
ison to data from [39J. The vertical line corresponds to the threshold for the 
string model (y/s = 2 GeV). Solid and dashed lines show results without and 
with qq annihilation channels (see subsection 12. 2|) . Errorbars on calculations are 
statistical. 
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Figure 2: Quark Diagram for the process ir + p — > Y, + K + . 
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Figure 3: Rapidity distributions of 7r ± and K ± for p+Be at 14.6 GeV/c in compar- 
ison to data from [T^j . Solid and dashed lines show results with energy dependent 
(see Eq.((TJ)) and constant (7 = 0.3) strangeness suppression factor, respectively. 
Errorbars on calculations are statistical. 
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Figure 4: Rapidity dependence of the inverse slope parameter T of the transverse 
mass spectra of ir + and K + for p+Be at 14.6 GeV/c in comparison to data from 
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Figure 5: Rapidity distributions of 7r ± and K for p+Au at 14.6 GeV/c in com- 
parison to data from jTJj- Solid and dashed lines show results with energy de- 
pendent (see Eq.((TJ)) and constant (7 = 0.3) strangeness suppression factor, re- 
spectively. The dotted line in the lower left panel shows a calculation without 
the qq annihilation channel in meson-baryon collisions. Errorbars on calculations 
are statistical. 
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Figure 6: Rapidity dependence of the inverse slope parameter T of the transverse 
mass spectra of tt + and K + for p+Au at 14.6 GeV/c in comparison to data from 
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Figure 7: Rapidity distributions of 7r ± and K ± for central collisions Si+Au at 14.6 
A GeV/c in comparison to data from pQ. Solid and dashed lines show results with 
energy dependent (see Eq.((T])) and constant (7 = 0.3) strangeness suppression 
factor, respectively. The spectra are divided by 28 in order to be able to compare 
to proton induced reactions. Errorbars on calculations are statistical. 
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Figure 8: Rapidity dependence of the inverse slope parameter T of the transverse 
mass spectra of 7r + and K + for Si+Au at 14.6 A GeV/c in comparison to data 
from pp. 
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Figure 9: Proton rapidity spectra from BUU events Au+Au at 10.7 A GeV double 
selected by the total multiplicity and the zero degree energy in comparison to data 
from Pj. The centrality decreases from the upper left corner to the lower right. 
The solid line shows the standard BUU calculation, whereas the dashed line is 
calculated with the medium modification described in section 
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Figure 10: Rapidity spectra of ir + from BUU events Au+Au at 10.7 A GeV 
selected by the zero degree energy in comparison to data from [9 . The centrality 
decreases from the upper left corner to the lower right. The solid line shows 
the standard BUU calculation, whereas the dashed line is calculated with the 
medium modification described in section 03 
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Figure 11: Rapidity spectra of K + from BUU events Au+Au at 10.7 A GeV 
selected by the zero degree energy in comparison to data from [9 . The centrality 
decreases from the upper left corner to the lower right. The solid line shows 
the standard BUU calculation, whereas the dashed line is calculated with the 
medium modification described in sectional 
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Figure 12: Transverse mass spectra for 7r + from Au+Au collisions at 10.7 A GeV 
for different slices of rapidity from the most central events selected by E Z cal- 
The rapidity slices range from 0.6-0.7 for the uppermost line to 1.3-1.4 for the 
lowermost line with a step 0.1. The spectra are multiplied by powers of 10: 
10°, 10 _1 , 10~ 7 from the uppermost to the lowermost line. Data are from 
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Figure 13: Transverse mass spectra for K + from Au+Au collisions at 10.7 A GeV 
for different slices of rapidity from the most central events selected by E Z cal- 
The rapidity slices range from 0.5-0.6 for the uppermost line to 1.2-1.3 for the 
lowermost line. The spectra are multiplied by powers of 10: 10°, 10" 1 , 10~ 7 
from the uppermost to the lowermost line. Data are from 9j. 
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Figure 14: Fiducial yields of tt + and K + divided by the number of projectile 
participants N pp as functions of N pp from Au+Au collisions at 10.7 A GeV in 
comparison to data from jH] . 
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Figure 17: Midrapidity yields for K as a function of energy in comparison to 
results of HSD, results of UrQMD and data from [201 IH1 12H] ■ 




Figure 18: Midrapidity yields for A + S as a function of energy in comparison 
to results of HSD, results of UrQMD and data from [HI EJ E2J EHJ E] . 
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Figure 19: K + /tt + ratio at midrapidity as a function of the beam energy. Upper 
left panel - comparison of the BUU results with the results of the UrQMD and 
HSD model |]|. Upper right panel - comparison between BUU and statistical 
model [22] • Lower left panel - comparison of standard BUU (cascade mode) 
and BUU with nuclear mean field potential. Lower right panel - comparison of 
the standard BUU calculation with the calculations without qq annihilation and 
without meson-meson collisions. The data are from [3J EJJ |2B] • 
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Figure 20: The K~ / ' K + ratio at midrapidity as a function of energy in comparison 
to results of HSD, results of UrQMD and data from [23 El 1221- The errorband 
on the BUU results indicate the uncertainty in the meson-meson channels (see 
discussions in sections ED and EJ). 
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Figure 21: The (A + S°)/7r ratio at midrapidity , where 7r = 1.5(7r + + 7r~), as a 
function of energy in comparison to data from j^OJ El H21 HH IS] • 
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Figure 22: Inverse slope parameter T for K + as a function of energy in comparison 
to data from [HI E31 • 
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Figure 23: Contribution of different channels to strangeness production in 
Au+Au collisions with b = 1 fm at 4 A GeV (upper panel), 10.7 A GeV (middle 
panel) and 20 A GeV (lower panel). 
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Figure 24: Time evolution of the central baryon density in Au+Au collisions at 
10.7 A GeV, b=0 fm. 
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Figure 25: In-medium modification factor F of the NN cross section vs cm. 
energy. Upper panel shows F at p = po for various outgoing channels: NNn - 
solid line, NNtttc - long-dashed line, NNp - short-dashed line, NAK - dotted 
line and NNKK - dash-dotted line. On the lower panel the factor F is presented 
for the NNn channel at p = p , 2p and 3p in the order from the upper to the 
lower line. 
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